ON THE MODELS OF NONLOCAL NONLINEAR OPTICS. 



BORIS G. KONOPELCHENKO AND ANTONIO MORO 

Dipartimento di Fisica dell'Universita di Lecce 
and Istituto Nazionale di Fisica Nucleare, Sezione di Lecce 
via Arnesano, 1-73100, Lecce 
E-mail: konopel@le.mfn.it, antonio.moro@le.infn.it 



We show that under certain assumptions a general model of nonlocal nonlinear 
response in 1 + 1— dimension is equivalent to the model considered by Krolikowski 
and Bang for a Kcrr-type medium. We derive the limit of weak nonlocality in high 
frequency regime and discuss the integrable cases. 



1. Introduction 

A paraxial laser beam propagating in a medium with response of the form 
D = Eq'E + cr'^D^^^ (cr << 1) is described by the following equation 

2zV^c.?^ + ViE + c.2D(3) =0, (1) 
oz 

where D'-'^^ is a certain function of the electric field. For derivation of such 
type of equation see [1, 2]. Recent theoretical and experimental results 
highlight a very interesting phenomenology, such as "accessible solitons" , 
collapse arrest, optical vortices (see e.g. [3—6] and references therein), asso- 
ciated with so-called nonlocal nonlinear media. Most of the models studied 
before are not integrable and exact analytical solutions are not available 
for them. However, recently, it was shown that certain choice of nonlocal 
responses leads to 2 -t- 1-dimensional integrable models 

In the present paper, for sake of simplicity we discuss the 1 + 
1— dimensional case. In the section 2, we consider a more fundamental 
nonlocal nonlinear model than one considered in the paper [7] . We specify 
under which conditions it leads us to a generalization of the model discussed 
by Krolikowski and Bang ^. The model proposed in [9] is of interest since 



•This work is partially supported by COFIN PRIN SINTESI 2004. 



2 



it can be solved analytically for weak nonlocality and it admits soliton-like 
solutions. 

In the section 3 we analyze a class of weak nonlocal responses in the 
high frequency regime. The study of high frequency limit can help to detect 

interesting properties of the beam such as singular phases and consequent 
vortex type behaviours. Moreover, this study could be useful to construct 
new nontrivial ansatz for applying a variational method. Looking for so- 
lutions which depend "slowly" on the coordinate along the propagating 
direction, we show that the phase of the electric field is given by an over- 
determined system of partial differential equations for the phase. Their 
compatibility and some integrable cases are also discussed. 



2. The nonlocal nonlinear Schrodinger equation. 

Equation (1) in 1 + 1— dimensional case looks like as follows 

We set £0 = 1 without lost of generality. Let us assume a constitutive 
relation of the following general form 

r + ca 

D(3) = / R{x-x';a) N{I{x'}) E{x^)dx'. (3) 



Distribution R{x — x':a) characterizes the nonlocal response around the 
point x and a is the "width" parameter (in the following it will be assumed 
to be depending on the frequency w). N{I) is an arbitrary nonlinear re- 
sponse depending on the intensity of the electric field / = |E|^. In what 
follows we will propose a way to simplify the general model (3) which will 
lead us to model considered in the paper [9] . 

Let us consider a nonlocal distribution R{x — x';a) of width 6R defined 
as the minimum such that 

R{x-x']a):^0, Vx' ^ [x - SR,x + 6R] . (4) 

Analogously wc; introduce the widths SE and SN of the electric field and 
the nonlinear response respectively. Suppose they verify the following con- 
ditions 

6R ~ 5N, 5R« SE. (5) 
Let us assume a nonlinear response of the form 

N{I{x))=NiX), (6) 
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where X = ^ x and 7 = l/SN. Expanding E(a;') and N {I{x')) in Taylor 
series around x, one gets the following approximation of the formula (3) 

D(3) ~ (^J^°^ R{x- x'; a) N {I{x')) dx'^ E(a;), (7) 

where we kept into account that due to the equation (6) higher orders of 
the expansion of N {I{x)) are not negligible. Note that the formula (7), 
in the case of nonlocal Kcrr-typc medium, leads to the nonlocal nonlinear 
Schrodinger equation discussed in the paper [9]. 

For instance, given a bell-shape electric field E = Eq exp [— a;^/2cr^], a 
nonlinear response of the form N{I) = 7« = |Eor" exp \- {-ixy /2(j'^\, 

where 7 = -\/2a, satisfies the condition (6). Nevertheless, it's easy to see 
that validity of relations (5) is sufficient to obtain the model (7). For 
instance, choosing E = Eo/cosh^(x) and (/) = condition (5) is 
verified for a large enough. 



3. Weak nonlocality and high frequency limit. 

Let us consider a limit of small nonlocality {5R « 1) such that it is 
reasonable to expand both E and N in Taylor series. This expansion gives 

B=NiIix))Eix) + R.[N'^ + 'J^E) + 

^ /^l^^^E IdNd^'E Id^Ndll ld^N\ 
^ \6 dx-^ 2 dx dx"^ 2 dx"^ dx 6 dx^ ) 

where the distribution R{x — x'\a) is assumed to be normalized 

+00 

R{x- x'] a) dx' = 1 (9) 



and Rn are the n-moments 

r'+oo 



/+00 
R{x-x';a) {x-x'f dx' . (10) 
-00 

In the present section, wc arc interested to perform the high frequency limit. 
In particular, we assume that for uj^ 00 the effective nonlocality decreases 
because of rapid oscillation of the electric field. So, one has 

lim R{x — x' ■,a{ijj)) = 5 {x — x') , (11) 
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where the parameter a is a function of the frequency such that there exists 
its hmit for w ^ 00 and 5{x — x') is the Dirac ^-function. 

As iUustrative example, we focus on the distribution defined as follows 



R{x — x';a) = < 



{-a^{x - x')'^ + a) x' e 



1 1 

-,x + 



a 



a 







x' ^ 



1 1 

X ,x-\ 

a a 



(12) 



It is straightforward to verify that distribution (12) tends to the J— function 
as a ^ oo, that is 



/oo 
R{x — x'; a) dx' = 1 
-oo 

/ + 00 
R{x-x';a) f{x') dx' = f{x). 
-oo 



A direct calculus shows that 



Rin 



r2n 
n2n' 



R2n+i =0, n e N U {0} . 



(13) 



(14) 



We would like to stress that the present discussion still holds for any non- 
local response whose moments are of the form (14). 

Let us consider a general dispersion law of the following form 



N = No + 



N 



(15) 



1 + {iroLo) 

where tq is usually called relaxation time. In the high frequency limit 
a; — > oo, expanding I = Iq + oj~'^'^Ii + . . . , one gets 

1 



N = No{Io) + ^'^y'K o 



,2u 



(16) 



where iVi is complex-valued and 2v > 0. For 21^ < 1, relation (15) is known 
as Cole-Cole dispersion law^° which substitutes the "classical" Debye law 
{21/ = 1) for a wide class of liquid and solid polar media. 

We perform the high frequency limit in usual way looking for solutions 
of the form 



E = Eo e' 



iujS 



(17) 



Moreover, we assume that the electric field depends slowly on the 
.2— variable according to the rule 



d_ 
dz 



'A 
dz 



(18) 
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It is straightforward to sec that under the assumption given above and 
under the choice a = 1 +v, the high frequency limit of equation (2) leads, 
at u"^ and w^"^"^ orders, to the following system 

Sl=No 

ro 1 (19) 

s, = -!1noSI + In„ 

where fx = df /dx. 

System (19) is over-determined and its compatibility condition is equivalent 
to the following nonlinear equation 

±N^^ No, = -r^No Nox + N^.. (20) 

For Ni^O solutions of equation (20) are complex- valued, and they implies 
complex- valued phases. Such solutions are connected in principle with ab- 
sorption or amplification of the electric field. 

We note that, once R is assigned, the model (7), at the order u)'^~'^^ , 
provides us with the equation 

= ^iVi, (21) 

and the compatibility condition can be written down in the form of conser- 
vation law Mqz = Nix where A^o = ±\/M). 

If we assume A^i = 0, setting Nq = Lp^ equation (20) is reduced to the 
well known Burgers-Hopf equation 

Vz = ±V> V>x- (22) 

It is solvable with the hodograph method and the solutions are given in 
terms of the following implicit relation 

x±ipz + H{ip)=Q, (23) 

where H is an arbitrary function of its argument. 

Now, let us assume A''i = N\ {Nq). Even in this case, equation (20) is 
reduced to the following quasilinear equation 

ATo, = ±N^ {N[ (No) - r2No) Nox, (24) 

where N[ = dNi/dNo. 

Analogously to the Burgers-Hopf equation it is solved by the hodograph 
relation 

X ± N^ {Ni (No) - r2No) z + L {No) = 0, (25) 
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where L is an arbitrary function of iVg. 

Other integrable equations can be obtained considering nonlocal re- 
sponse of the following form 

R{x-x';a)= i ^-{-^r + '^n^ ^ ' (26) 

i 0, x' 

where n e N and D = [x — x + i] . In this case, the n— moments are 

For any even n one gets the system (19). If n is odd, setting a= 1 + 2i/ we 
obtain the system 

Sl=No 

i n 1 (28) 

S. = -^N,S.,+ '-N,. 

Compatibility condition is equivalent to the equation 

Noz = i^r-iiVoiVoc. ± Vn^Ni^. (29) 
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